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Isospin effect of the statistical emission fragments from the equilibrated source is investigated in 
the frame of statistical binary decay implemented into GEMINI code, isoscaling behavior is observed 
and the dependences of isoscaling parameters a and [3 on emission fragment size, source size, source 
isospin asymmetry and excitation energies are studied. Results show that a and /3 neither depends 
on light fragment size nor on source size. A good linear dependence of a and 13 on the inverse 
of temperature T is manifested and the relationship of a = AC sym[{Z a / As)\ — [Zs I As)^lT and 



(i = ACsy^\(NslAs 



{Na/Aa)2\/T from different isospin asymmetry sources are satisfied. The 
symmetry energy coefficient Csym extracted from simulation results is ~ 23 MeV which includes both 
the volume and surface term contributions, of which the surface effect seems to play a significant 
role in the symmetry energy. 

PACS numbers; 24.10.-i, 24.10.Pa, 25.70.-z 
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I. INTRODUCTION 

The growing interest in isospin effects in nuclear reac- 
tions is motivated by an increasing awareness of the im- 
portance of the symmetry term in the nuclear equation 
of state. The availability of beams with large neutron-to- 
proton ratio, N/Z, provides the opportunity to explore 
the symmetry energy in very isospin-asymmetric nuclear 
systems. In such reactions, isospin degree of freedom has 
prominent roles and can be served as a valuable probe 
of the symmetry energy term of the nuclear equation of 
state. The isotopic composition of the nuclear reaction 
products contains important information on the role of 
the isospin on the reaction process. N/Z degree of free- 
dom and its equilibration, as well as the isospin asym- 
metry dependent terms of the nuclear equation of state 
(EOS) [ll,[3,[l,i,[Ell, have motivated detailed measure- 
ments of the isotopic distributions of reaction products. 

One important observable in heavy-ion collisions for 
determining the symmetry energy experimentally is the 
fragment isotopic composition investigated with the re- 
cently developed isoscaling approach [Tj, |8| . The isoscal- 
ing approach attempts to isolate the effects of the nu- 
clear symmetry energy in the fragment yields, thus al- 
lowing a direct study of the symmetry energy term in 
the nuclear binding energy during formation of hot frag- 
ments. Isoscaling refers to a general exponential relation 
between the yield ratios of given fragments between two 
reactions which differ only in their isospin asymmetry 
(N/Z). In particular, if two reactions, 1 and 2, lead to 
primary fragments having approximately same tempera- 
ture but different isospin asymmetry, the ratio R2i{N, Z) 
of the yields of a given fragment [N, Z) from these pri- 
mary fragments exhibits an exponential dependence on 
the neutron number N and the atomic number Z by fol- 



lowing form: 



R2{N, Z) ^ 



Y2{N, Z) 
Yi{N, Z) 



= Cexp(aAr-f /?Z), (1) 
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where a and (3 are two scaling parameters and C is 
an overall normalization constant. This scaling behav- 
ior has been observed in a very broad range of reactions 
[i,[i3,[n|Jl|Jl|J3jl^^ theoretical calcu- 

lationsHPlla^IIlll^illilli^ A 

review paper can be found in recent international "World 
Consensus Initiative" book on "Dynamics and Thermo- 
dynamics with Nuclear Degrees of Freedom" [s^] ■ 

The aim of the present paper is to investigate the 
isospin effect in sequential binary decay implemented into 
GEMINI code by isolating formation of the excited com- 
posite system, to learn the sensitivity of the isoscaling 
parameters with respect to the observable characterizing 
the source's state by only observing the first step decayed 
fragments, such as N/Z distribution and isotopic char- 
acteristics of the evaporated lightparticle. Sequential 
binary decay model GEMINI [sl, 0, [11] has been suc- 
cessfully used to describe the light particle evaporation, 
complex fragment emission, and N/Z distribution of the 
equilibrated compound source. 

The article is organized as follows. Section |TT] makes 
a simple review on the sequential binary decay imple- 
mented into GEMINI code, where a brief description 
of the symmetry energy adopted in the binding energy 
calculation and the method to extract temperature are 
given. In Section [1111 isoscaling phenomenon and its de- 
pendence on compound source size, excitation energy and 
isospin asymmetry are presented, the sensitivities of the 
isoscaling parameters a and (3 to observables character- 
izing the state of the source are discussed in detail. The 
contribution from surface effect on the symmetry energy 
term is discussed in Section IIVI Finally a summary is 
given. 



2 



II. MODEL OVERVIEW 

GEMINI model [H, Hi] calculates the decay of com- 
pound nuclei by modes of sequential binary decays. All 
possible binary divisions from light-particle emission to 
symmetric division are considered. The model employs a 
Monte Carlo technique to follow the decay chains of indi- 
vidual compound nuclei through sequential binary decays 
until the resulting products are unable to undergo further 
decay. 

The decay width for the evaporation of fragments with 
Z < 2 is calculated using the Hauser-Feshbach formalism 
[3^ . For the emission of a light particle (Zi,Ai) with 
spin Ji from a system {Zo,Ao) with excitation energy 
E* and spin Jg, leaving the residual system ^2,^2 with 
spin J2, the decay width is given by: 



T{Zi,Ai,Z2,A2) = 



2Ji + l 
27rpo 



J0+J2 

E 

= |./0-./2 



L 



E'-B-Erot{J2) 



Ti{e)p2{U2,J2)de. 



(2) 



In this equation / and e are the orbital angular momen- 
tum and kinetic energy of the emitted particle, P2(C^2, J2) 
is the level density of the residual system with thermal 
excitation energy 

U2=E* -B~Erot{J2)-£, (3) 

where B is the binding energy, Erot{J2) is the rotation 
plus deformation energy of the residual system, po is level 
density of the initial system and T; is the transmission 
coefficients. 

For binary divisions corresponding to the emission of 
heavier fragment, the decay width is calculated using the 
transition state formalism of Moretto (37| 

V{Zi,Ai,Z2,A2)^-^^^ / Psad{Usad,Jo)de, 



Jo 



(4) 

where Usad and psad are the thermal energy and level 
density of the conditional saddle-point configuration, re- 
spectively, 



Usad — E* ~ Esad{Jo) — £, 



(5) 



where Egad is the deformation plus rotation energy of 
the saddle-point configuration and e now is the kinetic 
energy of the translational degree of freedom. 

The symmetry energy term due to the neutron-proton 
excess is presented in calculating the masses of nuclei. 
For heavy systems (Z>12), the masses of the initial 
and residual systems are obtained from the Yukawa-plus- 
exponential model of Krappe, Nix and Sierk [s^ without 
the shell correction, pairing correction term for odd-odd 
nuclei is included. The parameters for this model are 
taken from the fit to experimental masses of MoUer and 
Nix [3^. For very light systems (A<12), masses of the 
nuclei are calculated from the experimental ones. 

Taking into account the effect of a predicted increase 
in the symmetry energy associated with the temperature 
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FIG. 1: (Color online) Isoscaling parameters a (positive val- 
ues) and f3 (negative values) as a function of the fragment 
proton number Z or neutron number N from source pair {Zs 
= 50, Ab = 100) and {Zs = 50, As = 105) at various excita- 
tion energies: Eex = 1.0 (open squares), 1.4 (solid squares), 
2.0 (open circles), 2.4 (solid circles), 3.0 (open up triangles) 
MeV/nucleon. 

dependence of effective nucleon mass in the surface of the 
nucleus [4l|, the kinetic part of the symmetry energy is 
related to the level spacing at the Fermi surface and so 
it is also related to the level density parameter ot, the 
following temperature-dependent kinetic symmetry en- 
ergy was therefore included to calculate the temperature 
dependent level density parameter ut in the GEMINI 
simulations (ssl. [ioj 

i?,^™(r) = 0.82247(— --)(iV-Z)2, (6) 
" ar ao 

where ao = (1.6^ + 1.8A^/^)/15.5. 

The effective thermal excitation energy of the equi- 
librated system before light particle evaporation and 
saddle-point configuration before heavy fragment emis- 
sion can be got by 



E. 



the 



= U2-E: 



kin 
sym 1 



or 



ther 



Usad ~ E' 



kin 



(7) 



(8) 



Then the nuclear system temperature is approximately 
Eiher/a^. (9) 

III. ISOSCALING BEHAVIORS 

To make a systematic study of source parameters 
which might influence the isoscaling behavior, in our 
present work several pairs of equilibrated sources are con- 
sidered at various initial excitation energies Ef.x — 1-0, 
1.4, 2.0, 2.4 and 3.0 MeV/nucleon. To avoid possible ef- 
fects of different magnitudes of Coulomb interaction on 
isotopic distributions, we consider pairs of sources with 
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FIG. 2: (Color online) Isoscaling parameters a (positive val- 
ues) and 13 (negative values) as a function of the fragment 
proton number Z or neutron number A*' from source pairs 
with the fixed proton number Zs —50 at excitation ener- 
gies Eex = 2 MeV/nucleon. Symbols in figure correspond 
to Ya^^ii5/Ya,=iio (solid squares), Ya,=iio/Ya,=io5 (open 
square), yA^=io5/'K4,=ioo (solid circles), Ya,=ii5/Ya,^io5 
(open circles), Ya^=iio /Ya, =ioo (solid up triangles). 
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FIG. 3: (Color online) Same as Fig. 2 but for source 
pairs with the fixed proton number Zs =30 at excitation 
energies Eex ~ 2 MeV/nucleon. Symbols in figure cor- 
respond to YA^=69/YA,^6(i (solid squares), YA^^6(i/YA^=63 
(open square), 1a, =637X4, =60 (solid circles), Ya, =69/1^^=63 
(open circles), Ya^^66/Ya^=60 (solid up triangles). 

the same proton numbers Zs but different mass numbers 
As- The equilibrated source pairs are chosen at differ- 
ent mass region and system isospin asymmetry N/Z, two 
groups of the source pairs have been used: (1) Zg = 50 
with As = 100, 105, 110 and 115, respectively; (2) Zs = 
30 with As — 60, 63, 66 and 69, respectively. We adopt 
the widely used convention to denote with the index "2" 
the more neutron-rich system and with the index "1" the 
more neutron-poor system. In this situation the value of 
a is always positive because more neutron-rich clusters 
will be produced by the neutron-richer source and the 



value of /? is always negative. The yield ratios R2i{N, Z) 
are calculated and the corresponding isoscaling behaviors 
are investigated over all possible decayed fragments. In 
the present study, we will focus on the first decay step 
to investigate the decay products even though a real de- 
cay chain is usually longer than one for the deexcitation 
process of a thermal source. The advantage of studies 
on the first step is that the source size and temperature 
parameters are well defined. 

Fig. [1] plots the isoscaling parameters a and /3 as a 
function of fragment proton number Z and neutron num- 
ber N for source pairs with fixed proton number: {Zs 
= 50, As = 100 {Ns/Zs = 1.0) and As = 105 (Ns/Zs 
= 1.1)) at excitation energies 1.0, 1.4, 2.0, 2.4 and 3.0 
MeV/nucleon. Fig. [5] and Fig. [3] plot the isoscal- 
ing parameters a and /3 as a function of fragment pro- 
ton number Z and neutron number N for different size 
source pairs at excitation energies 2.0 MeV/nucelon re- 
spectively. Fig. [2] is the source pairs with source atomic 
number Zs — 50 and different mass number As — 100, 
105, 110 and 115, and Fig. [3] is for the source pairs with 
source atomic number Zs = 30 and different mass num- 
ber As = 60, 63, 66 and 69. 

Fig. [TJ [5] and [3] show that the values of isoscaling 
parameters a and /3 are essentially flat with the frag- 
ment proton number Z or neutron number N in the case 
ploted, as well as cases which included in the simulations 
which were not plotted in the figure samples. Average 
isoscaling parameters a and P can be calculated over the 
flat region to discuss the dependence of a and /3 on the 
properties of emission source, such as excitation energy, 
source size and source asymmetry of the isospin. The 
average a is calculated over the range Z < 9 and the 
average /3 is calculated over the range < 9, which 
keep both average values are calculated over on the same 
fragments region. In the following sections isoscaling pa- 
rameters a and /3 refer to the calculated average values. 



A. Dependence of a and f3 on excitation energy 

The excitation energy and temperature dependences of 
a and /3 are shown in Fig. 3] and Fig. [5] for different size 
and isospin asymmetry source pairs, respectively. Tem- 
perature of the initial source are tabulated in TABLE [H 
in which the level density parameter ax is derived from 
the Eq. ^ and the temperature are calculated from Eq. 

Both parameters a and \f3\ in Fig. [Jand Fig. [5] have 
monotonic dependence on the excitation energy and their 
absolute values decrease with the excitation energy, \(3\ 
generally are greater than a. In right panel of Fig. 3] 
and Fig. [5l a and \f3\ show linear dependence on 1/T, 
which evidences the relationship of a = A/i„/T and 
(3 = A^p/T [2l|. The slope of the relation between a 
(1/3 1) and 1/T should be the free neutron (proton) chem- 
ical potential difference A/i„ (|A/ip|), the linear depen- 
dence of a (|/3|) on 1/T also evidences the constant of 
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FIG. 4: (Color online) Dependence of isoscaling parame- 
ters Q and |/3| on excitation energy (left panel) and the 
inverse temperature 1/T (right panel) for various source 
pairs with same proton number Zg = 30. Symbols in the 
figure are a (solid symbols) or (open symbols) from 
source pair Ya^=63/Ya^^60 (squares), Ya^^66/Ya^=63 (cir- 
cles), YAs=69/YA^^a6 (up-triangles), Va^ =66/^^^=60 (down- 
triangles), and Yas=69/Yas^63 (diamonds). 
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FIG. 5: (Color online) Same as Fig. 4 but for vari- 
ous source pairs with ~ 50. Symbols in the figure 
are a (solid symbols) or \f3\ (open symbols)from source 
pair 1a,=io5/'K4.,=ioo (squares), Ya^=iio/Ya,=io5 (circles), 
Ya,=ii5/Ya^^iio (up-triangles), yA,=iio/yA,=ioo (down- 
triangles), and Yas=ii5 /Ya^^io5 (diamonds). 

free neutron (proton) chemical potential difference A/i„ 
(|A/j,p|) between two initial sources with defined asym- 
metry isospin, which are independent of the excitation 
energy or temperature. The linear dependence of isoscal- 
ing parameters a and /3 with 1/T offers a signal to cal- 
culate the free neutron (proton) chemical difference A/i„ 
(|A/ip|), which can be calculated from the slope fitting in 
the right panels of Fig. 2] and Fig. [51 Considering a > 
and /3 < 0, the free neutron (proton) chemical potential 



TABLE I: Level density parameters ar and temperature T 
for different source systems and excitation energies. 



Source 




(MeV) ar 


TiMeVy- 


Zs^SO As = 


60/63 


1.0 


7.5/7.8 


2.8/2.9 






1.4 


7.3/7.6 


3.4/3.5 






2.0 


7.1/7.4 


4.1/4.2 






2.4 


7.0/7.3 


4.5/4.6 






3.0 


6.8/7.1 


5.1/5.2 


Zs^30,As= 


66/69 


1.0 


8.1/8.5 


2.9/2.9 






1.4 


8.0/8.3 


3.5/3.5 






2.0 


7.7/8.0 


4.2/4.2 






2.4 


7.6/7.8 


4.6/4.7 






3.0 


7.4/7.6 


5.2/5.3 


Zs^50,As= 


100/105 


1.0 


11 8/19 3 


2.8/2.9 






1.4 


11.5/11.9 


3.4/3.5 






2.0 


11.0/11.4 


4.2/4.3 






2.4 


10.8/11.2 


4.6/4.7 






3.0 


10.5/10.9 


5.3/5.4 


Zs^50,As= 


110/115 


1.0 


12.8/13.2 


2.9/3.0 






1.4 


12.4/12.8 


3.5/3.6 






2.0 


11.9/12.3 


4.3/4.3 






2.4 


11.6/12.0 


4.8/4.8 






3.0 


11.3/11.8 


5.4/5.4 


"refers to the saddle-point configuration temperal 


ure 


TABLE II: 


Difference of free neutron and proton chemical 


potential fj,„i — fipi for series isospin asymmetry systems from 


the slope fitting in Fig. |4]and[5l 




System 






n - Mp'(MeV) 


Zs As^ As2 










30 60 63 


2.79±0.21 




3.42±0.33 


6.21±0.39 


30 63 66 


2.94±0.29 




3.53±0.30 


6.47±0.42 


30 66 69 


1.64±0.24 




2.42±0.53 


4.06±0.58 


30 60 66 


5.38±0.23 




5.77±0.49 


11.25±0.54 


30 63 69 


4.35±0.20 




4.71±0.22 


9.06±0.30 


50 100 105 


3.66±0.28 




2.96±0.29 


6.62±0.40 


50 105 110 


2.58±0.35 




2.40±0.10 


4.98±0.36 


50 110 115 


1.89±0.15 




2.13±0.29 


4.02±0.33 


50 100 110 


5.53±0.48 




5.41±0.49 


10.94±0.69 


50 105 115 


3.95±0.11 




4.06±0.53 


8.01±0.54 



"refers to fin^A,^ - Mn.A^j 
''refers to fip^A^^ - t^P.A,-^ 

"lefeis to (m„,a.,2 - /^p,As2 ) - (Mn.Aei - t^-p.A^-^ ) 



difference A^„ > and A^p < 0. 

For an isospin asymmetry system / = [N — Z)/A, 
its neutron and proton chemical potential are fj,n,i and 
fipj, respectively. From the isoscaling calculation we can 
get the free neutron and proton chemical potential dif- 
ference A^n = Unj^ - Unj^ and A^p = Hpj^ - Hpj^ , 

thus A^n - ^I^J-n = (Ain,/2 " A*n,/i ) " (Mp,/2 " /^P,/l)- ^ 

one isospin symmetry system with / = is included, 
fin, 1=0 = fJ'p,i=o, the difference of free neutron and pro- 
ton chemical potential /i„7 — /ip/ of the isospin asymme- 
try system can be deduced as shown in [ij, [l^. This 
^J'nI — pi-pi extracted from Fig. 2] and are printed in 
TABLE ini in which the index of isospin asymmetry pa- 
rameter / was replaced by the mass number Ag. 
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FIG. 6: (Color online) q-T (positive parts) and fi-T (negative 
parts) as a function of {ZslAa)\ — {Zs/As)2 or {Ns/As)i — 
{Ns/As)2 of the sources for various source pairs with atomic 
number Zs = 50 which is illustrated in the figure at different 
excitation energies. 
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FIG. 7: (Color online) Same as Fig. 6 but for various source 
pairs with Zs = 30 which is illustrated in the figure at different 
excitation energies. 

B. Dependences of a. and (i on isospin asymmetry 
and source size 

To explore the origin of isoscaling behavior and the de- 
pendence of isoscaling parameters a and /3 on the system 
size and the isospin composition, we performed calcula- 
tions on source systems with different sizes and asymme- 
try {Ns/Zs) values. Symmetry energy coefficient Csym is 
dependent on not only the nuclear density p, but also the 
system temperature [4J|. In GEMINI investigation, the 
nuclear density is set to be around the saturation density 
Pq so that the system density influences on the extraction 
of symmetry energy coefficient can be neglected in this 
case. The symmetry energy coefficient Csym can be seen 



only temperature dependent. Since the systems have dif- 
ferent sizes, the temperatures of the system are slightly 
different even though the excitation energy is fixed at 
same values. In this case, in Fig. [S]and[7]the temperature 
is served as a correction factor, we made the temperature 
correction by a • T and (i ■ T. Fig. [5] and [7] depicts a • T as 
a function of {Zs/As)\ — {Zs/As)\ and /3 • T as a function 
of {Ns/As)\ — {Ns/As)1) of the initial source pair, plot- 
ted in the same figure for different excitation energies. 
The source pairs we simulated are Zs — 50 pairs {Zs = 
50: As = 100, 105, 110 and 115) and Zs = 30 pairs {Zs 
= 30: As — 60, 63, 66 and 69). All these systems with 
different source sizes and isospin asymmetries lie along 
one line, which illustrates that the isoscaling parameters 
a and (3 are not sensitive to the system size and charge. 
The linear fits of the calculation points in Fig. [5] and [7] 
are printed on the figure, slopes in Fig. [6] are 92.0 ± 2.8 
and 95.0 ± 1.6 in Fig. [3 respectively. This approximate 
linear relationship has been also observed in many exper- 
imental data [l^, [ij] as well as other model calculations 
with EES, SMM model [IH, AMD [H and IQMD model 

M- 

It has been shown in the framework of the grand- 
canonical limit of the statistical multifragmentation 
model [2l| and in the expanding-emitting source model 
[l3| , a simple relationship can be derived between isoscal- 
ing parameters a and A{Zs/As)'^ , (3 a,ndA{Ns/Asy of 
the sources by following forms: 



, C sym 

a — 4:- 



and 



T 



^Csym 



-(Z^Y _ (Zs_ 

XaJi \a. 



Zs_Y 

2. 



T 



(10) 



(11) 



Above equations have been proved to be good approx- 
imations by many calculations and experimental data, 
and are generally adopted to constrain the symmetry en- 
ergy coefficient Csym in experiment. Here they are ver- 
ified by the results in Fig. [H] and [3 i.e. the slope of 
a ■ T with respect to A{Zs/As)'^ and f3 ■ T with respect 
to A{Ns/AsY. By finear fitting the data in Fig. [6]and[7] 
with Eq. pU]) and pT|) the symmetry energy coefficient 



Csym can be constrained to be 23.0 ± O.TMeV from Fig. 
H and 23.8 ± OAMeV from Fig. [71 comparable with the 
values from other models and experiments [l^.lT3l. [l5l . [2ll |. 



IV. 



SURFACE CONTRIBUTION TO ISOSPIN 
EFFECT 



GEMINI code involves the isospin effect or symmetry 
energy in calculating the binding energy. For heavy sys- 
tems {Zs >12), the mass excess of the initial and residual 
systems of a spherical nucleus is given by following equa- 
tion 
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FIG. 8: (Color online) Top panel: mass excess calculated 
from Eq. (I12|l as a function of nuclear mass number, curves 
from left to right correspond to the atomic number Z of 
isotopes series from Z = 5 to Z = 60 in order; Bottom 
panel: numerical calculated derivative of symmetry energy 
E sym (MeV/nucleon) with respect to neutron number A'^ in 
Eq. (|12p as a function of {Z/Af (solid marks) and the linear 
fit of the numerical calculation, fitted result is printed on the 
Figure. 



where / = {N — Z) /A. The parameters are taken from 
[11], M„ = 8.071431, Mp = 7.289034, = 1.4399764, b 
= 0.99, W = 36, aei 1.433e-■^ tq = 1.16, a = 0.68, a, 
= 21.13, ks = 2.3, ay = 15.9937, k^ = 1.927, A = l2/^/A, 
5 — 20/ A. Correction arising from single-particle effect 
is added to Eq. (fT2)) . which includes Coulomb energy for 
diffusive surface, proton form factor and charge asymme- 
try term [H, [sl. El]. In the top panel of Fig. [51 the 
calculated mass excess by Eq. ([T^ is plotted as a func- 
tion of mass number A, the isotopes with different charge 
number Z from 5 to 60 are presented by lines from left 
to right. 

In Eq. p2p the dominant symmetry energy term 
arises from two parts: the volume symmetry energy part 
ttykyl^A and the surface diffuseness term related with the 
isospin asymmetry —askgl^ {A^^^ + ...... | . As derived in 

other models and theoretical frame [1, |2l[ , isotope yield 
ratio is dominantly determined by the symmetry term 



in the binding energy for two equilibrium sources with 
comparable mass and temperature but different isospin 
degree Ng/Zg. In this case, the isoscaling parameter a 
can be achieved by following approximate form 

a = -As„/T, (13) 

where As„ is the difference in neutron separation energy 
between the two sources, considering the dominant term 
in separation energy is symmetry term, which is calcu- 
lated in GEMINI simulation by Eq. (fT3)) . The symmetry 
term taken from Eq. can be expressed alone by 



Eg 



Cyl^A - 



Cgl'{ 



A 



2/3 



(14) 



the first term in Eq. is the volume term of isospin 
asymmetry part, the second term is the surface effect 
of isospin asymmetry. Cy = Qyky = 15.9937 x 1.927 — 
30.8199 MeV which is the generally used symmetry 
energy coefficient at saturate nuclear density po=0.16 
/m-3, cs = asks = 21.13 x 2.3 = 48.599 MeV. The differ- 
ence in neutron separation energy between two sources 
can be approximately obtained by taking the deriva- 
tives of the symmetry energy of Eq. with re- 
spect to N. Numerical calculation of the derivatives of 
Eq. HI]) with respect to TV (i.e. AEsym/ AN) is per- 
formed, and approximate linear function on (Z/A)'^ is 
observed, shown in the bottom panel of Fig. [51 The lin- 
ear fits give AEsym/AN = 22.26019 - 90.5205(Z/A)2 w 
22.63(1 — 4(Z/A)^), effective symmetry energy coefficient 
Csyrn in the equation of state deduced here approximates 
22.63 MeV, which is consistent with the symmetry energy 
coefficient Csym derived from isoscaling parameters a and 
j3 very well, namely Csym = 23.0±0.7 MeV derived from 
a or Csym = 23.8 ± 0.4 MeV from /3. As we know, the 
derivatives with respect to TV of the volume term in Eq. 
(0 is c„(l - 4(Z/v4)2) = 30.8199 - 123.2796(Z/A)2 = 
30.8199(1 - 4(Z/y4)2), hence the contribution from the 
surface effect is approximately 22.630125(1 - /^(ZjAf) - 
30.8199(1 -4(Z/A)2) = -8.55971(1-4(^/^)2). gum of 
the volume and surface terms of the symmetry energy 
give s the consistent result with the experimental results 
[2l| and other model calculations, which indicates that 
the surface effect term of the symmetry energy strongly 
influences the isoscaling parameter a and /3. 
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V. SUMMARY 

The statistical sequential decay process of the equili- 
brated source has been successfully investigated by the 
GEMINI model. The model constrains the source size 
and fragments density at saturate nuclear density po and 
the decay starts from an equilibrated system which is 
separated from any dynamical process. By only calcu- 
lating the first step decay, the source size, density and 
temperature can be well defined. In our work we apply 
it to survey the evolution of isospin degree of freedom. 
Isospin effects on the decayed fragments from different 
source sizes, isospin asymmetries and excitation energies 
have been systematically investigated. Isoscaling phe- 
nomena are observed for the emission fragments of light 
fragments. Isoscaling parameters a and |/3| decrease with 
the increasing of excitation energy, they generally show 
linear dependence on the inverse of the system temper- 
ature T, and are independent of the source size. From 
the linear function between a versus 1/T or /3 versus 
1/T, the extracted linear slope is the free neutron (pro- 
ton) chemical potential difference A/i„ and A/ip between 
two systems. If the isospin symmetry source is included, 
then the difference of neutron and proton chemical po- 



tential A/i„ — A/ip can be determined, thus the slope 
difference between a and /? versus 1/T offer a signal to 
calculate the neutron and proton chemical potential dif- 
ference. After considering system temperature, the sim- 
ple relationship a = \C sym{{ZY I A^f - {Z2/A2f]/T and 
= '^Csym[{N-i_/Axf - lN2/A2f]/T can be well repro- 
duced for various sources. Only when the surface effect 
term in the symmetry energy is taken into accounted, the 
above linear relationship between a and A(Z/A)^, or (3 
and A(A^/^)^ can give consistent symmetry energy co- 
efficient Cgym with experimentally proposed results and 
other model results, which illustrates that the surface ef- 
fect plays a significant role in the symmetry energy term. 
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